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1. Introduction 

One of the outstanding problems in science today is the identity of the dark matter of the uni- 
verse The existence of dark matter is implied by a large number of observations, including 
the dynamics of galaxy clusters, the rotation curves of individual galaxies, the abundances of light 
elements, gravitational lensing, and the anisotropics of the cosmic microwave background radia- 
tion. The energy density fraction of the universe in dark matter is 23%. The dark matter must be 
non-baryonic, cold and coUisionless. Cold means that the primordial velocity dispersion of the dark 
matter particles is sufficiently small, less than about 10^^ c today, so that it may be set equal to zero 
as far as the formation of large scale structure and galactic halos is concerned. CoUisionless means 
that the dark matter particles have, in first approximation, only gravitational interactions. Particles 
with the required properties are referred to as 'cold dark matter' (CDM). The leading CDM candi- 
dates are weakly interacting massive particles (WIMPs) with mass in the 100 GeV range, axions 
with mass in the 10^^ eV range, and sterile neutrinos with mass in the keV range. 

Today I argue that the dark matter is axions ^. The argument has three parts. First, axions 
behave differently from the other forms of cold dark matter because they form a Bose-Einstein 
condensate [Q]. Second, there is a tool to distinguish axion BEC from the other forms of CDM 
on the basis of observation, namely the study of the inner caustics of galactic halos. Third, the 
evidence for caustic rings of dark matter is consistent in every aspect with axion BEC, but not with 
WIMPs or sterile neutrinos. 

Before I start, let me mention that H. Baer and his collaborators have shown that in many 
supersymmetric extensions of the Standard Model, the dark matter is axions, entirely or in part [Q]. 



2. Axions 

Shortly after the Standard Model of elementary particles was established, the axion was pos- 
tulated to explain why the strong interactions conserve the discrete symmetries P and CP. For 
our purposes the action density for the axion field (p{x) may be taken to be 

where m is the axion mass. The self-coupling strength is 

in terms of the axion decay constant / and the masses /m„ and rud of the up and down quarks. In 
Eq. (^n|), the dots represent higher order axion self-interactions and interactions of the axion with 
other particles. All axion couplings and the axion mass 



10^2 GeV 

T 



m ~ 6 • 10"^ eV " (2.3) 



are inversely proportional to /. / was first thought to be of order the electroweak scale, but its 
value is in fact arbitrary [^. However, the combined limits from unsuccessful searches in particle 
and nuclear physics experiments and from stellar evolution require / > 3 • 10^ GeV |^. 



2 



The dark matter is mostly an axion EEC 



Pierre SIKIVIE 



Furthermore, an upper limit / < 10^^ GeV is provided by cosmology because light axions are 
abundantly produced during the QCD phase transition In spite of their very small mass, these 
axions are a form of cold dark matter. Indeed, their average momentum at the QCD epoch is not 
of order the temperature (GeV) but of order the Hubble expansion rate (3 • 10^^ eV) then. In case 
inflation occurs after the Peccei-Quinn phase transition their average momentum is even smaller 
because the axion field gets homogenized during inflation. For a detailed discussion see ref. [^. In 
addition to this cold axion population, there is a thermal axion population with average momentum 
of order the temperature. 

The non-perturbative QCD effects that give the axion its mass turn on at a temperature of order 
1 GeV. The critical time, defined by m{ti)ti = 1, is fi ~ 2 • 10"^ sec (//lO'^ GeV) 3. Cold axions 
are the quanta of oscillation of the axion field that result from the turn on of the axion mass. They 
have number density 

4-10^^ / / \'(aihy ' 
cm^ V 10^^ GeV J \ a{t) 

where a{t) is the cosmological scale factor. Because the axion momenta are of order ^ at time t\ 
and vary with time as fl:(f)^\ the velocity dispersion of cold axions is 

mt\ a[t) 

if each axion remains in whatever state it is in, i.e. if axion interactions are negligible. Let us 
refer to this case as the limit of decoupled cold axions. If decoupled, the average state occupation 
number of cold axions is 

(21)' { f 



Clearly, the effective temperature of cold axions is much smaller than the critical temperature 



for Bose-Einstein condensation. Bose-Einstein (BEC) may be briefly described as follows: if 
identical bosonic particles are highly condensed in phase space, if their total number is conserved 
and if they thermalize, most of them go to the lowest energy available state. The condensing 
particles do so because, by yielding their energy to the remaining non-condensed particles, the 
total entropy is increased. 



Eqs. ( |2.6[ ) and ( |2.7| ) tell us that the first condition is overwhelmingly satisfied. The second 
condition is also satisfied because all axion number violating processes, such as their decay to two 
photons, occur on time scales vastly longer than the age of the universe. The only condition for 
axion BEC that is not manifestly satisfied is thermal equilibrium. Thermal equilibrium of axions 
may seem unlikely because the axion is very weakly coupled. However, it was found in ref. [Q] 
that dark matter axions do form a BEC, marginally because of their self-interactions, but certainly 
as a result of their gravitational interactions. No special assumptions are required. 
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3. Bose-Einstein condensation of cold dark matter axions 

Axions are in thermal equilibrium if their relaxation rate F is large compared to the Hubble 
expansion rate H{t) = ^. At low phase space densities, the relaxation rate is of order the particle 
interaction rate = naSv where (J is the scattering cross-section. The cross-section for <p + <p — > 
(p + (p scattering due to axion self interaction is in vacuum 



1 A 



2 



m 



If one substitutes Gq for a, is found much smaller than the Hubble rate, by many orders of 
magnitude. However, in the cold axion fluid background, the scattering rate is enhanced by the 
average quantum state occupation number of both final state axions, a ~ Gq^^, because energy 
conservation forces the final state axions to be in highly occupied states if the initial axions are in 
highly occupied states. In that case, the relaxation rate is multiplied by one factor of ^ [^] 

r ~ « ao 5v ^ . (3.2) 



Combining Eqs. ( |2.4| - |2.qj3.l| ), one finds r{ti)/H{ti) ~ ^(1), suggesting that cold axions thermal- 
ize at time t[ through their self interactions, but only barely so. 

A critical aspect of axion BEC phenomenology is whether the BEC continues to thermalize 
after it has formed. Axion BEC means that (almost) all axions go to one state. However, only if 
the BEC continually rethermalizes does the axion state track the lowest energy state. 



The particle kinetic equations that yield Eq. ( p.2\ ) are valid only when the energy dispersion 
^m{dv)^ is larger than the thermalization rate [|l^]. After ti this condition is no longer satisfied. 
One enters then a regime where the relaxation rate due to self interactions is of order [Q, |n|] 

r;t ~ A « m^^ . (3.3) 

Tx{t)/H{t) is of order one at time t\ but decreases as t a{t)^^ afterwards. Hence, self interactions 
are insufficient to cause axion BEC to rethermalize after t\ even if they cause axion BEC at ti. 
However gravitational interactions, which are long range, do the job later on. The relaxation rate 
due to gravitational interactions is of order [|, 11] 



Y^^Gnn?f (3.4) 

where / ~ [mdvY^ is the correlation length. Tg{t) /H{t) is of order 4 • 10"^(//10'^ GeV)3 at time 
t\ but grows as ta^^{t) oc a{t). Thus gravitational interactions cause the axions to thermalize and 
form a BEC when the photon temperature is of order 100 eV (// 10^^ GeV) 2 . 

The process of axion Bose-Einstein condensation is constrained by causality. We expect over- 
lapping condensate patches with typical size of order the horizon. As time goes on, say from t to2t, 
the axions in f-size condensate patches rethermalize into 2?-size patches. The correlation length is 



then of order the horizon at all times, implying 5v ~ instead of Eq. (2.5), and Tg/H oc t a (t) 



after the BEC has formed. Therefore gravitational interactions rethermalize the axion BEC on ever 
shorter time scales compared to the age of the universe. The question now is whether axion BEC 
has impUcations for observation. 
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4. Dark matter caustics 



The study of the inner caustics of galactic halos [12, 13] may provide a useful tool. An isolated 
galaxy like our own accretes the dark matter particles surrounding it. Cold coUisionless particles 
falling in and out of a gravitational potential well necessarily form an inner caustic, i.e. a surface of 
high density, which may be thought of as the envelope of the particle trajectories near their closest 
approach to the center. The density diverges at caustics in the limit where the velocity dispersion of 
the dark matter particles vanishes. Because the accreted dark matter falls in and out of the galactic 
gravitational potential well many times, there is a set of inner caustics. In addition, there is a set 
of outer caustics, one for each outflow as it reaches its maximum radius before falling back in. We 
will be concerned here with the catastrophe structure and spatial distribution of the inner caustics 
of isolated disk galaxies. 

The catastrophe structure of the inner caustics depends mainly on the angular momentum 



distribution of the infalling particles ||13|]. There are two contrasting cases to consider. In the first 
case, the angular momentum distribution is characterized by 'net overall rotation'; in the second 
case, by irrotational flow. The archetypical example of net overall rotation is instantaneous rigid 
rotation on the turnaround sphere. The turnaround sphere is defined as the locus of particles which 
have zero radial velocity with respect to the galactic center for the first time, their outward Hubble 
flow having just been arrested by the gravitational pull of the galaxy. The present turnaround radius 
of the Milky Way is of order 2 Mpc. Net overall rotation implies that the velocity field has a curl, 
V X V 7^ 0. The corresponding inner caustic is a closed tube whose cross-section is a section of 



the elliptic umbilic (D_4) catastrophe [12, 13]. It is often referred to as a 'caustic ring', or 'tricusp 
ring' in reference to its shape. In the case of irrotational flow, V x v = 0, the inner caustic has a 
tent-like structure quite distinct from a caustic ring. Both types of inner caustic are described in 



detail in ref.0130. 

If a galactic halo has net overall rotation and its time evolution is self-similar, the radii of its 
caustic rings are predicted in terms of a single parameter, called jmax- Self-similarity means that the 
entire phase space structure of the halo is time independent except for a rescaling of all distances 
by R{t), all velocities by R{t)/t and all densities by l/t^ [0 [Bl, [T|, [T^. For definiteness, R{t) 



will be taken to be the turnaround radius at time t. If the initial overdensity around which the halo 
forms has a power law profile 

'-^^(±r , (41) 

where M,- and 5M,- are respectively the mass and excess mass within an initial radius r,-, then R{t)oc 



2,_2_ , , 

fs + gc [14|. In an average sense, e is related to the slope of the evolved power spectrum of density 



perturbations on galaxy scales [18]. The observed power spectrum implies that e is in the range 



0.25 to 0.35 [|l6|]. The prediction for the caustic ring radii is {n = 1, 2, 3, .. ) [ p^ , 17] 

^ 40kpC / V„t \ / jmax \ 

n V220km/sy' \0A%) 

where Viot is the galactic rotation velocity. Eq.( ^) is for e = 0.3. The a,, have a small e depen- 
dence. However, the a„ oc l/n approximate behavior holds for all e in the range 0.25 and 0.35, so 
that a change in e is equivalent to a change in jmax- (£i7max) = (0.30, 0.180) implies very nearly 
the same radii as (£,7max) = (0.25, 0.185) and (0.35, 0.177). 
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Observational evidence for caustic rings with the radii predicted by Eq. (4.2) was found in the 



statistical distribution of bumps in a set of 32 extended and well-measured galactic rotation curves 



[]19[], the distribution of bumps in the rotation curve of the Milky Way |gOp, the appearance of a 
triangular feature in the IRAS map of the Milky Way in the precise direction tangent to the nearest 
caustic ring [|20|], and the existence of a ring of stars at the location of the second (n = 2) caustic ring 
in the Milky Way [^T]]. Each galaxy may have its own value of jmax- However, the j'max distribution 
over the galaxies involved in the aforementioned evidence is found to be peaked at 0.18. There is 
evidence also for a caustic ring of dark matter in a galaxy cluster [^]. 

Recently the rotation curve of our nearest large neighbour, the Andromeda galaxy, was mea- 



sured with far greater precision and detail than hitherto achieved [gSp. The new rotation curve has 
three prominent bumps, at 10 kpc, 15 kpc and 29 kpc. The positions of these bumps are in the 
ratios predicted by the caustic ring model and thus provide fresh additional evidence. 

5. The caustic ring halo model 

The caustic ring model of galactic halos [|l^] is the phase space structure that follows from 
self-similarity, axial symmetry, and net overall rotation. Self-similarity requires that the time- 
dependence of the specific angular momentum distribution on the turnaround sphere be given by 

Em 

i{n,t)=j{n)^ (5.1) 

where n is the unit vector pointing to a position on the turnaround sphere, and j{n) is a dimen- 
sionless time-independent angular momentum distribution. In case of instantaneous rigid rotation, 
which is the simplest form of net overall rotation, 

/(«) = jmax nx{zxn) (5.2) 



where £ is the axis of rotation and jmax is the parameter that appears in Eq. (4.2). The angular 
velocity is w = iasi.£ Each property of the assumed angular momentum distribution maps onto an 
observable property of the inner caustics: net overall rotation causes the inner caustics to be rings. 



the value of jmax determines their overall size, and the time dependence given in Eq. ( |5.1| ) causes 
a„ oc l/n. 

The angular momentum distribution assumed by the caustic ring halo model may seem implau- 
sible because it is highly organized in both time and space. Numerical simulations [ ^4] ] suggest that 
galactic halo formation is a far more chaotic process. However, since the model is motivated by 
observation, it is appropriate to ask whether it is consistent with the expected behaviour of some 
or any of the dark matter candidates. In addressing this question we make the usual assumption, 
commonly referred to as 'tidal torque theory', that the angular momentum of a galaxy is due to 
the tidal torque applied to it by nearby protogalaxies early on when density perturbations are still 



small and protogalaxies close to one another [25, 26]. We divide the question into three parts: 1. is 



the value of jmax consistent with the magnitude of angular momentum expected from tidal torque 
theory? 2. is it possible for tidal torque theory to produce net overall rotation? 3. does the axis of 



rotation remain fixed in time, and is Eq. (5.1) expected as an outcome of tidal torque theory? 
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6. Magnitude of angular momentum 

The amount of angular momentum acquired by a galaxy through tidal torquing can be reliably 
estimated by numerical simulation because it does not depend on any small feature of the initial 
mass configuration, so that the resolution of present simulations is not an issue in this case. The 
dimensionless angular momentum parameter 

X = ^\ , (6.1) 

where G is Newton's gravitational constant, L is the angular momentum of the galaxy, M its mass 
and E its net mechanical (kinetic plus gravitational potential) energy, was found to have median 



value 0.05 [[27|]. In the caustic ring model the magnitude of angular momentum is given by jmax- 
As mentioned, the evidence for caustic rings implies that the j^ax -distribution is peaked at j^^x — 
0.18. Is the value of jmax implied by the evidence for caustic rings compatible with the value of A 
predicted by tidal torque theory? 

The relationship between jmax and A may be easily derived. Self-similarity implies that the 



halo mass M{t) within the turnaround radius R{t) grows as t^'e p4|]. Hence the total angular mo- 
mentum grows according to 

dt J dQdt 9e 

where we assumed, for the sake of definiteness, that the infall is isotropic and that j{n) is given by 



Eq. (5.2). Integrating Eq. (6.2), we find 



Similarly, the total mechanical energy is 

J R{t) dt 5-3e R{t) 

Here we use the fact that each particle on the turnaround sphere has potential energy —GM{t) /R{t) 



and approximately zero kinetic energy. Combining Eqs. (|6.1|), ( |6.3| ) and ( |6.4| ) and using the relation 
R(tf = ^fGM{t) g, we find 



For £ = 0.25, 0.30 and 0.35, Eq. (Q imphes A / jmax = 0.281, 0.283 and 0.284 respectively. Hence 
there is excellent agreement between jmax — 0.18 and A ^ 0.05. 

The agreement between jmax and A gives further credence to the caustic ring model. Indeed 
if the evidence for caustic rings were incorrectly interpreted, there would be no reason for it to 
produce a value of jmax consistent with A. Note that the agreement is excellent only in Concordance 
Cosmology. In a flat matter dominated universe, the value of jmax implied by the evidence for 



caustic rings is 0.27 [12, 17|. 
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7. Net overall rotation 



Next we ask whether net overall rotation is an expected outcome of tidal torquing. The answer 
is clearly no if the dark matter is coUisionless. Indeed, the velocity field of coUisionless dark matter 
satisfies 

^(r,0 = |^(?,0 + (v(r,0-V)v(r,0 = -V(/.(r,0 (7.1) 
where <p{r,t) is the gravitational potential. The initial velocity field is irrotational because the 



expansion of the universe caused all rotational modes to decay away [28]. Furthermore, it is easy 
to show | [T3| ] that if V x v = initially, then Eq. ( [7?T| ) implies V x v = at all later times. Since net 
overall rotation requires V x v 7^ 0, it is inconsistent with coUisionless dark matter, such as WIMPs 
or sterile neutrinos. If WIMPs or sterile neutrinos are the dark matter, the evidence for caustic 
rings, including the agreement between j^iax and A obtained above, is purely fortuitous. 

Axions [^, 1^ ^ differ from WIMPs and sterile neutrinos. Axions are not coUisionless, in 
the sense of Eq. (fO|), because they form a rethermalizing Bose-Einstein condensate. This process 



is quantum mechanical in an essential way and not described by Eq. ([TJJ). By rethermalizing we 
mean that thermahzation rate remains larger than the Hubble rate so that the axion state tracks the 
lowest energy available state. The compressional (scalar) modes of the axion field are unstable 
and grow as for ordinary CDM, except on length scales too small to be of observational interest 
[0]. Unlike ordinary CDM, however, the rotational (vector) modes of the axion field exchange 
angular momentum by gravitational interaction. Most axions condense into the state of lowest 
energy consistent with the total angular momentum, say L = Lz, acquired by tidal torquing at a 
given time. To find this state we may use the WKB approximation because the angular momentum 
quantum numbers are very large, of order 10^'' for a typical galaxy. The WKB approximation maps 
the axion wavefunction onto a flow of classical particles with the same energy and momentum 
densities. It is easy to show that for given total angular momentum the lowest energy is achieved 



when the angular motion is rigid rotation. So we find Eq. ( p. 21 ) to be a prediction of tidal torque 
theory if the dark matter is axions. 

Thermahzation by gravitational interactions is only effective between modes of very low rela- 
tive momentum because only in this case is the correlation length /, that appears in Eq. (p.4[), large. 



After the axions fall into the gravitational potential well of the galaxy, they form multiple streams 



and caustics like ordinary CDM [29]. The momenta of particles in different streams are too differ- 
ent from each other for thermahzation by gravitational interactions to occur across streams. The 
wavefunction of the axions inside the turnaround sphere is mapped by the WKB approximation 
onto the flow of classical particles with the same initial conditions on that sphere. The phase space 
structure thus formed has caustic rings since the axions reach the turnaround sphere with net overall 
rotation. The axion wavefunction vanishes on an array of lines. These lines, numbering of order 
10^*', may be thought of as the vortices characteristic of a BEC with angular momentum. However, 
the transverse size of the axion vortices is of order the inverse momentum associated with the radial 
motion in the halo, (mv^)^' ~ 20 meters for a typical value (10^^ eV) of the axion mass. In a BEC 
without radial motion the size of vortices is of order the healing length [^], which is much larger 
than (mv,-)^^ 
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One might ask whether there is a way in which net overall rotation may be obtained other than 
by BEC of the dark matter particles. I could not find any. General relativistic effects may produce a 
curl in the velocity field but are only of order (v/c)^ ~ 10^^ which is far too small for the purposes 
described here. One may propose that the dark matter particles be coUisionfuU in the sense of 
having a sizable cross-section for elastic scattering off each other. The particles then share angular 
momentum by particle collisions after they have fallen into the galactic gravitational potential well. 
However, the collisions fuzz up the phase space structure that we are trying to account for. The 
angular momentum is only fully shared among the halo particles after the flows and caustics of 
the model are fully destroyed. Axions appear singled out in their ability to produce the net overall 
rotation implied by the evidence for caustic rings of dark matter. 

If the dark matter is WIMPs or sterile neutrinos, the velocity field of dark matter is curl-free. 
As already mentioned, the inner caustics of galactic halos have then a tent-like structure which is 



quite distinct from caustic rings [|13|]. Also, the angular momentum of the dark matter accreting 
onto a halo is not shared among the infalling particles. The total angular momentum vector L of 
the halo is the same as for axion dark matter, since it is determined by the outcome of tidal torque 
theory, but unlike the axion case is the sum of many contributions randomly oriented with respect 
to one another. The tent-like inner caustics have therefore random orientations, whereas the caustic 
rings of the axion case he all in the galactic plane. 

8. Self-similarity 

The third question provides a test of the conclusions reached so far. If galaxies acquire their 
angular momentum by tidal torquing and if the dark matter particles are axions in a rethermalizing 
Bose-Einstein condensate, then the time dependence of the specific angular momentum distribution 



on the turnaround sphere is predicted. Is it consistent with Eq. (p.l|)? In particular, is the axis of 
rotation constant in time? 

Consider a comoving sphere of radius S{t) = Sa{t) centered on the protogalaxy. a{t) is the 
cosmological scale factor. S is taken to be of order but smaller than half the distance to the nearest 
protogalaxy of comparable size, say one third of that distance. The total torque applied to the 
volume V of the sphere is 

Z{t)= [ d^r5p{r,t)rx{-V<p{r,t)) (8.1) 
Jv{t) 

where dp{f,t) = p{f,t) — po{t) is the density perturbation. po{t) is the unperturbed density. In 
the linear regime of evolution of density perturbations, the gravitational potential does not depend 
on time when expressed in terms of comoving coordinates, i.e. (p{r = a{t)x,t) = <p{x). Moreover 
5{r,t) = has the form 5{r = a{t)x,t) = a{t)5{x). Hence 

T{t) = pQ{t)a{tf I d^x 5{x)xx (-V,0(x)) . (8.2) 
Jv 



Eq. (8.2) shows that the direction of the torque is time independent. Hence the rotation axis is time 
independent, as in the caustic ring model. Furthermore, since Po{t) a{t)^ , z{t) oc a{t) oc ?i and 
hence £{t) oc L{t) oc Since R{t) oc ^a+w, tidal torque theory predicts the time dependence of 
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Eq. (5.1) provided e = 0.33. This value of e is in the range, 0.25 < £ < 0.35, predicted by the 
evolved spectrum of density perturbatuions and supported by the evidence for caustic rings. So the 
time dependence of the angular momentum distribution on the turnaround sphere is also consistent 
with the caustic ring model. 



9. Conclusion 



If the dark matter is axions, the phase space structure of galactic halos predicted by tidal torque 
theory is precisely, and in all respects, that of the caustic ring model proposed earlier on the basis of 
observations. The other dark matter candidates predict a different phase space structure for galactic 
halos. Although the QCD axion is best motivated, a broader class of axion-Uke particles behaves 
in the manner described here. 
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